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Abstract 



In the adiabatic approximation, most of the effects of qq loops on spec- 
troscopy can be absorbed into a static interquark potential. I develop a 
formalism which can be used to treat the residual nonadiabatic effects as- 
sociated with the presence of nearby hadronic thresholds for heavy quarks. 
I then define a potential which includes additional high energy corrections 
to the adiabatic limit which would be present for finite quark masses. This 
"improved" potential allows a systematic low energy expansion of the impact 
of thresholds on hadronic spectra. 
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I. INTRODUCTION 



The valence quark model is surprisingly successful at describing mesons and baryons 
as qq and qqq systems moving in effective potentials. The surprise comes in part because 
hadrons are so strongly coupled to their (real and virtual) decay channels that each nearby 
channel ought to shift a hadron's mass by Am ~ F typical, thereby totally disrupting the 
valence quark model's spectroscopy. 

A simple resolution of this conundrum has been proposed in a series of papers |l],|2j 
examining the effects of "unquenching the quark model", i.e., allowing extra qq pairs to 
bubble up in valence quark states. This bubbling dresses the valence hadrons with a certain 
class of meson loop diagrams f3[. (These papers also address how the OZI rule |4| survives 
unquenching; in this paper I will exclusively consider flavor nonsinglet states for which such 
OZI- violation is not an issue.) The proposed resolution is an extension of the idea [||] that 
in the absence of light quarks the heavy quarkonium potential y^ diabatlc ^ r s j ^ b Qr j s the 
adiabatically evolving ground state energy E Q (r) of the purely gluonic QCD Hamiltonian in 
the presence of a static color triplet source Q and color anti-triplet sink Q separated by a 
distance r . Once rif light quarks are introduced into this Hamiltonian, two major changes 
occur: 

1. E (r) will be shifted to E nf (r) by ordinary second order perturbation theory, and 

2. E Hf (r) will no longer be isolated from all other adiabatic surfaces: once pair creation 
can occur, the QQ flux tube can break to create states {Qq) a {qQ)p with adiabatic 
energy surfaces that are constant in r at the values e a + ep (e^ is the i th eigenvalue of 
the Qq system, with the heavy quark mass tuq subtracted). 

Despite the latter complication, in the weak pair creation limit the flux-tube-like adiabatic 
surface E nf {r) can be tracked through the level crossings that occur when E n Ar) = e a + ep 
and identified as the renormalized QQ adiabatic potential V^f mbaUc (r) = E nf (r). In Ref. 
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it is shown that for large r, y^ ahatlc (y^ remains linear, so that the net effect of the pairs 
is simply to renormalize the string tension. Since quark modellers determined their string 
tension from experiment, the quark model heavy quarkonium potential already included the 
effect of meson loops to leading order in the adiabatic approximation, i.e., b n . = b, the 
physical string tension. 

Note that a similar renormalization occurs at short distances: in lowest order — > 
af 1 ^ = 127r/[(33 — 2nf)ln(Q 2 / Aq CD )}. The renormalization of the string tension by qq loops 
is quite similar, though complicated by the existence of the open channels corresponding to 
adiabatic level crossings. It should also be stressed that the possibility of subsuming qq 
loops into b n . only occurs if one sums over a huge set of hadronic loop diagrams (real and 
virtual) Jl| . No simple truncation of the sum over loops, as is often attempted in hadronic 
effective theories, is generally possible. Consider, for example, the simplest orbital splitting 
a2(1320) — p(770). Summing the mass shifts associated with the known decay modes of these 
states would significantly change their absolute masses and violently alter their splitting. 
Preserving them requires a large renormalization of the string tension and summing over 
loop graphs involving many high mass (i.e., virtual) channels, since qq creation inside the 
original QQ state is dual to a very large tower of (Qq) a (qQ)p intermediate states. 

Although the renormalization y^ diabattc _ » yadiabatic ca pture the bulk of the effect of 
"unquenching" in heavy quarkonia, E n (r) deviates quite substantially from linearity near 
level crossings |l| . Both this fact and explicit modelling suggest that for phenomenologically 
relevant quark masses substantial nonadiabatic effects will remain after renormalization, 
and in particular that states near thresholds to which they are strongly coupled should be 
expected to deviate from their potential model positions. This paper is devoted to developing 
a method for addressing these residual effects. This is straightforward as txiq — > 00, but I 
will show that for finite mq it is essential to go beyond the naive adiabatic approximation 
to define an "improved" interquark potential which includes the high energy part of the 
corrections to the adiabatic limit. 
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II. THE FORMALISM IN THE ADIABATIC LIMIT 



To deal with violations of the adiabatic approximation, we can closely imitate the normal 
methods of mass renormalization. For very massive quarks Q and Q, the effects of all 
hadronic loop graphs can be subsumed into 

y adiabatic ^y,^ yadiabatic _|_ ^ ' ^y adiabatic 

where v^ odm6o * JC (r) is the "purely gluonic" static QQ potential, and ^y^g iabatlc (r) is the shift 
in this static potential generated by the channel a(3. Here the subscript on Vq atlc i s used 
to denote that it is purely gluonic; we have suppressed additional labels to identify which 
gluonic adiabatic surface y^ diabatlc represents (the normal meson surface, the first A = ±1 
hybrid surface, etc.) since our discussion applies identically to them all. For the low-lying 
thresholds of interest to us here, /\y^ iabatlc ^ typically have a strength of order Aqcd 
and a range of order Aqq D . This range arises because ^aifqo) an d ^p{f q Q) are localized 
at relatively small |r*qQ | and |r„<g| for low- lying states so that for large r the production of 
such states by the point-like creation of a qq pair is strongly damped by the rapidly falling 
tails of their confined wavefunctions; conversely, for small r the created q and q are easily 
accommodated into the "heart" of their respective wavefunctions. 

Let us now compare the adiabatic Hamiltonian for the QQ system [[J 

TT P I adiabatic /n\ 

-[^adiabatic 7> ' nt \ I 

(with the reduced mass of rrii and rrij) with the two channel Hamiltonian H^ a ^ that is 
the penultimate step in generating H a( aabatic in the sense that all channels except a[3 have 
been integrated out: 



#(«/3) 



p I '[/'adiabatic A T/ adiabatic TJll 

2 Mqq + % V ^ n W) 



,,2 



H (*P) 2^ + € « + e P 



(3) 



where ffL is an interaction which couples the QQ system to the single channel (Qq) a (qQ)/3 
with the matrix elements dictated by the underlying pair creation Hamiltonian H^. In the 
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adiabatic limit we must recover H ac nabatic from H^>, but ifw 3 ) contains the full dynamics 
of the coupling of the QQ system to the channel a(3. With the superscript on H^ a ^ we are 
making explicit that if ( a/3 ) has the channel a(3 removed from the QQ adiabatic potential 
and added back in full via Hf^y We could in general remove any subset of n channels from 
y adiabatic ^ an( j ac ^ ^hem back in dynamically as part of an (n + 1)— channel problem. In 
the limit of taking all channels we would recover the original full unquenched Hamiltonian. 
However, since our treatment is in lowest-order perturbation theory, the effects of the indi- 
vidual channels are additive, and Eq. ([$[) with just an individual channel {af3) selected for 
study is sufficient for our purposes. 

Note that the hadronic multichannel version of our unquenched Hamiltonian is an ap- 
propriate representation of qq pair creation in a confined system. When the pair is created, 
the {QqqQ) system has three relative coordinates which we may take to be p, the separation 
between the center of mass of meson (3 and that of meson a, and the two intrameson coor- 
dinates rqQ = rq — rQ and r q Q = f q — fg. Since we ignore the residual final state interaction 
between the color singlets (Qq) a and (qQ)p, the eigenstates of this sector are mesons in 
relative plane waves, corresponding to the entry H%2 i n Eq. (§). Thus, with p p canonically 
conjugate to p, the three quantum labels (p p ,a,/3) replace the three labels (p,r q Q,r q Q). 

The main goal of this paper is to describe the relation between the eigenvalues of the 
adiabatic Hamiltonian (0) and the dynamic Hamiltonian (|3]). If we define 



H 

and 



p I if adiabatic n 
2MQQ ^ »/ 



+ e a + e/3 



(4) 



Hpert 



A T r adiabatic TJll 
-^ V a{3 n (af3) 

(a/3) U 



(5) 



and denote the QQ eigenvalues of iio and H^ a/3 ^ by Ef and E\ a ^\ respectively, then since 
fl(ap) _ jj^ _|_ ii pert) by second order perturbation theory AE^ 01 ^ = E^' — E® is given by 
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^ piadiabatic(al3) ^ ^j^dynamic(a/3) ^y-^ 

where I^q) is the i th eigenstate of H . This simple equation is the main focus of this paper. It 
represents the correction to the adiabatic approximation for the QQ energy eigenvalues from 
a full dynamical versus an adiabatic treatment of the channel (aft). In what follows I will 
first show explicitly that AE\ a ^ — * as expected in the limit mq — > oo. I will then define an 
improved effective potential v^ 1 J iproved which incorporates "trivial" high energy corrections 
to the adiabatic approximation, but which is essential for incorporating threshold effects in 
a systematic low energy expansion for finite mg. 

I begin by defining precisely /W^ tabaUc in Eq. (|]). If \aj3(p)) denotes an a/3 state with 
relative coordinate p, then as niq — > oo 

(aM\Rg\QQW) = (a(3(p)\HZ\QQ(r)) = c a p(f)5 3 (p - f) (8) 

since the QQ relative coordinate is frozen in the adiabatic approximation by definition and 
since tq — tq — > p as niQ — > oo. Thus 



br - (e Q + ep) 



J br-(e a + e l3 ) 
= 5 3 {f' - r)AV$ iabatic (r) (11) 
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for br far from e a + e^; see Ref. for a discussion of how the poles in (|T0| ) are to be handled. 

radi 
a/3 



Given this AV^ iabatic {f), by definition 



AE adiabaUc(af3) = ^^yadiabatic^ ^ 

. (13) 



/ a/3 

br - (e a + ep) 

I now show how AEf habaUc ^ a ^ approximates the true shift 

AE dynamic(aP) = j, 3q I (<*ff(g) I #gfl 1^) ? ^ 



Ef-(e a + ep + 



even for "nearby" thresholds as mg — > oo. Denote by (v)i the expectation value of the 
variable v in the state \i/)q). In the limit tuq — ► oo, each of 2 ^ ^ , 6(r)j, and vanishes 
like ( A ^ CP Y^J^qcd and so is small compared to e a + ep, which is of order Aqcd, but large 
compared to the corrections to e a + e^, which are of order Aqcd/^q- (In the general 
power law potential c n r n , they each behave like ( Aqcd ) n ^ n+2 AQ C D, i-s., they vanish for any 

2 

confining (n > 0) potential). For ; this statement is nontrivial: it relies on the behavior 
of the numerator of Eq. (0) . Using Eq. (§) , 

(a^lflJ^igQ®) = 7335/ dV e ^>%3(r) (15) 

= c a/3 (p- q) , (16) 

so, even though \p\ ~ (Ag^m^) 1 / 3 — > 00, |p — g| must be of order Aqcd since c a/ 3 is a light 
quark object. After writing E® = & ^ + 6(r)j, we can therefore Taylor series expand: 

dynamic(af3) , 1 s f J% „\ I of -*\\ Trqq |„/,i\l 2 



(1 + (— ^- + 6(0*] + • • •) (17) 

e« + C/3 •/ •/ •> 
(l + (^—)[^^ + ^] + ■•■), (18) 

where 4> l (p} = — / dPre'^ip^r). Noting that 



rf's c a/3 (s) = c a/3 (0) (19) 

and that except for the c the factors of the integrand are slowly varying functions, we can 
approximate 

~c a p(s) ~ c a/3 (0> 3 (s) + ■ • ■ (20) 

to obtain 

AE dynamic(aP) „ f ^^2^ + ^_J_ ^)i ~ + ^ r ).] + . . .) (21) 

e a + C/3 J £a + C/3 m Q 
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f d*M(P)\ 2 il + i^—)} (22) 



e a + eg J ' e a + 6p ' 

F br- (e a + e/9 ) 



&r-(e Q + e/3 ) 

This expression differs slightly from ^^ adm6atlf w) j n gq ^ has |c a/ g(r)| 2 — > |c Q( g(0)| 2 . 

However, for low-lying states, ( j^|gjj 2 )i — 1 ~ ( A ^ fl ) 2 ^ 3 which is negligible as m Q — > oo 
compared to ( - + )i ~ ( A ^ fl ) 1 '' 3 which we retained. (The physics behind this approxima- 
tion is simply that |c aj g(r)| 2 reflects light quark scales while iV'o^l 2 reflects short distance 
scales as tjiq — > oo.) thus to leading order as mg — > oo, Ai^ ' = AE i ' for 

low-lying states, as we set out to prove. 

The deviation AE\ a ^ of the energy of the state i from its value in the adiabatic potential 
yadiabatic ^ ue ^ Q ^ e residual dynamical effects of the channel (Qq) a (qQ)p thus has the 
property that AE\ a ^ — > channel by channel as itlq — > oo. The utility of this formalism 
for heavy quarkonium arises from not only this property, but also that it allows one to 
consistently focus on low-lying thresholds. The latter feature is based on the fact that the full 
shift AEi = J2 a p AE\ may receive significant "random" contributions from strategically 
placed low mass channels, but for fixed large itlq the AE^^ — ► rapidly as e a + gets 
large. This rapid convergence occurs because, first, as e a + e@ — ► oo, Ef and g 2 /2/i aj g in Eq. 
([yD and 6r in Eq. (|13|), which were already small with respect to e a + ep even for low mass 
channels, become negligible. In this limit AE\ a ^ is therefore trivially zero independent 
of the accuracy of the approximations inherent in Eqs. (|17[)-(P4|). Moreover, the factor 
|c Q( g(r)| 2 in the numerator of each of AE° ,diabatlc ^ a ^ anc [ ^^jd,ynamic(a/3) ra pj^y approaches 
zero as e a + ep — > oo since for low-lying states \ip®) there is little kinetic energy in the initial 
wavefunction and the pair creation process can only create momenta of order Aqcd- We 
can therefore expect this formalism to provide a rapidly converging low-energy expansion of 
the effects of pair creation on heavy quarkonia. 
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III. AN IMPROVED QUARKONIUM POTENTIAL 



While suitable for heavy quarkonia, the framework of Section II has serious shortcom- 
ings for light quark spectroscopy. Because the eigenvalues ei m< ^ and eg™ and the matrix 
elements (a^ ft m v\p)\Hf c \QQ{r)) for finite m<g are only qualitatively related to their 
rriQ — > oo counterparts e a , ep, and (a(3(p)\H^\QQ(f)) , Av^f abaiiC (r) is not in this case an 
accurate approximation to the effects of the channel {Qq) a {(lQ)f3- As a result, the AE\ a ^ 
will not be small, i.e., the critical separation of the effects of the channel aft into large 
adiabatic and small residual dynamical effects will fail. If it were only for a few nearby 
channels, this failure would not be so serious, but while Ef and g 2 /2/i Q/ g in the finite mg 
analogue of Eq. ([14] ) and br in Eq. ( |13D can still be neglected as e a + ep — > oo, since 
+ ef Q) + e a + e p and (a^ (3^\p)\Hf c \QQ(f)) ± (aft(p)\Hf c \QQ(f)}, the finite 

i £ a T-idynamicallaB) . nn , , . . n 1 a i-iadiabatic(al3) . , 1 • 1 • 1 

m,Q analogue of AE i will not trivially approach AE i m this nigh energy 

limit. Moreover, while the matrix elements (a^ mQ ^/3^ mci '(p)\H^\QQ(r)) may still be ex- 
pected to cut off high mass channels, since the momenta in low mass states and in the pair 
creation process are comparable, these channels will be cut off more slowly in the finite rriQ 

i £ a T-id,ynamical(aP) ,1 . a T-i&diabaticlotB) mi 1 j. 1 j_i 

analogue of AE i than m l\E i . lnese shortcomings make the vtlq — > oo 

2/3 ■ 



framework far less useful in light quark systems since AEi = J2 a s AE\ a ^ will converge only 



marginally faster than the "brute force" sum J2 a /3 /\Ef ynamc ^ al3 ' > . 

I will now show that it is possible to define an improved effective quarkonium potential 
yimproved ^jgjj l eac l s to energy shifts 5E\ a ^ which vanish as e a + ep — > oo for any itlq and so 
give smaller and more rapidly converging corrections to the quark model spectroscopy built 
on y^P TOVed than the AE\ al3 \ The price to be paid for this important feature is that the 
universal (flavor- independent) adiabatic quarkonium potential y^ iabatlc must be replaced by 
a flavor-dependent effective potential v^ iproved built out of y^ diabatlc pl us flavor-dependent 



j • l j • t rirnproved 

contributions Av a/3 



The basic idea is very simple. For any mq ||, the shift in the energy of the state l^o ^ mQ ^) 
due to channel a(3 is given by the generalization of Eq. fll^l), namely 
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. j^dynamic{aP) {m Q ) _ f .3 J (^g^ (g ggg) ^ ^) P 

where the superscripts (?71q) denote quantities at finite uiq in contrast to those previously 
defined for vtlq — ► oo. 

I begin by examining the limit ei™' 3 '* + e^ mQ ' ) >> e\ (mQ ' ) and 5^—, eacn °f which are 
for itiq comparable to Aqcd themselves of order Aqqd- I n this limit we have 

("»q) (mg) 

where 

-1 

( m Q) 



A^ ( ; Q) = {mQ) {mQ) J ^^)l« (mQ) ^ (WQ H^)(« (mQ) ^ (mQ H^I^&) (27) 
is an mg-dependent but \ip l ^^-independent effective potential operator. Thus 



-1 



ea v ' + e 



/3 



(a^/S^J^I^igQCf)) (28) 



d 3 q d 3 p ' <?p (QQ{r')\Hf a0) \a^(3^\p')) 



(m Q ) , (m Q ) / " * " ^ " ^ VW 



^-(a^^^lflSnlQgCTO) (29) 



6a + £g 



(« (mQ) /3 (mQ) (p)|^)|gO(r)) (30) 



which can be compared to Eq. (|9|). I next introduce the finite uiq analogue of Eq. (^j). 
As rug — ► oo, the form of Eq. (§) is model independent, with dynamical effects residing in 
c a/ g(r). However, for finite mg even the form of the analogue to Eq. (§) becomes model- 
dependent. The key to extending the utility of the heavy quarkonium framework down to 
light quark masses is to make a "local approximation" 
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( a MpM(fl\Hg\QQ(fi) * ^\r)5\p- rf) (31) 
where as mg — > oo, r\ — » 1 and 

^(0-^(0, (32) 



the right hand side being the function defined in the adiabatic limit by Eq. (U). Note 
that (r) involves at the microscopic level overlap integrals between \QQ(r)} and 
|a;( ,n o)/3( m «)(^)^ with wavefunctions vk?) and ^"^(^q) for finite mg, while c a p(r) 

involves the heavy quark limits ijJaifqo) an d ^(^q) °f these wave functions. In the sim- 
plest and most common models 0,0,0, the "local approximation" is exact and automatic 



with f] = mQ ® m , corresponding to qq pair creation that is point-like and instantaneous. 



There are, of course, other possibilities, both local and nonlocal [ IT5 1 ; the latter need to be 
"localized" via an approximation of the form of Eq. (|3TD in order that one may define their 
improved effective potentials. 

With © used in ©, we have 

(QQ(r')\AV^\QQ(r)) = 5*(r'-f>)\ g£ ] (33) 



e a + e 



= S 3 (f'- f)AV^ proved (r) . (34) 
By construction, ^y^-P roved gives a AEl mproved< " a ^ which converges to AEf ynamic ^ a/3 ^ mQ ^ for 



approximation (|31~D is made, ^y^P roved w jU gj ve a m0 re accurate approximation to the 
high energy effects of than AV^3 ia6atlc (f), so its use leads to a more rapidly converging 



models in which the local approximation ( pT|) is exact. In all cases |L0 |, once the local 

a/3 

approximation to the effects of thresholds on the hadronic system. We can actually improve 
matters even further if we incorporate the additional convergence available even for low mass 
thresholds as rriQ — ► oo by defining || 



rj 3 \c [ Z Q \r]r)\ : 



a -irimproved i -a _ / vr_n /qr\ 

^ = 7 / (m ) 7 (m Q )J • ^ 

or - (ea v + v ) 
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( m <j) i (me) 



For eT Q > + e p 

improved i 



large, br is negligible so this approximation is no worse than that leading 



to AV r 



(f), but ^v^ provea (r) also approaches our old AV^ abattc (r) as vtlq 
so gives a good approximation to AEf ynamtc ^ a/3 ^ mQ ^ f or a \\ a p j n ij m jt. 



oo and 



Given these features, we can improve upon Eqs. (H) and (§) by defining 

,2 



if, 



improved 



P I jrimproved 
2/iQQ 



(36) 



and 



(a/3) 



improved 



p^ |_ -^improtied _ improved 



2 ^QQ 



rrqq 



rrqq 



(37) 



along with the analogue of Eq. (^) 



a/3 



Wo 



) + ae; 



dynamic(a/3) (mq ) 



(38) 



The <^ a/3) now approach zero both in the strict adiabatic limit tuq — > oo for all a/5 and 



also in the limit e 



(m Q ) (m Q ) 



+ e 



/3 



>> Aqcd for all mg. They therefore allow a systematic low 



energy expansion of the impact of thresholds on the spectra of all quarkonia. 



IV. CONCLUSIONS 

I have presented here a formalism for calculating the nonadiabatic component AE" 13 of 
the mass shift of a valence heavy QQ state % from the hadronic loop process % — *> af3 — > i, 
i.e., the component of this process that cannot be absorbed into the renormalized heavy 
quarkonium potential. The resulting formula was shown to have the expected property that 
as rriQ — > oo. The formula is also very simple and, when combined with a pair 
creation model like the flux-tube-breaking model 0] or the 3 Po model ||, should provide 
a quick method of estimating the influence of nearby thresholds on the spectra of heavy 
quarkonia. 

I have also shown how to define an "improved" quarkonium potential which incorporates 
nonadiabatic effects associated with high mass thresholds for any vtlq. When this potential 
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is identified with the empirical quark model potential, the deviations 5E\ a ^ of the spectrum 
from the potential model predictions due to thresholds have the property that they vanish 
both as rriQ — > oo for all a[3 and also as the mass e^ Q ^ + e^ 1 ^ of the threshold a(3 gets large 
for any itlq. This improved potential therefore allows a systematic low energy expansion of 
the impact of thresholds on hadronic spectra. 

This advantage has a price: the "improved" potential has the characteristic that it 
violates the rule of flavor independence. While this rule is valid in the heavy quark limit 
and to leading order in perturbative QCD for light quarks, violations are to be expected. 
Indeed, though obscured by possible relativistic corrections, there are indications from quark 
models that the best empirical potentials are system-dependent ]ll . 



An important step not taken here is to calculate the AE^ al3 ^ and 8E\ a ^ for selected 
channels to assess numerically how rapidly each converges as eiT^ + tg 1 — > oo Q, and 
to quantify the mQ-dependence of v^" iproved . Quark models seem to constrain this mass 



dependence to be surprisingly weak [p^j . Assuming that the approach defined here passes 
quantitative tests such as these, it will then be interesting to apply it to a number of 
outstanding phenomenological issues. Among these are the threshold shifts in the cc and bb 
systems and the A(1520) — A(1405) problem. It will also be amusing to study heavy-light 
systems to see explicitly how groups of states conspire to maintain the spectroscopic relations 



required by heavy quark symmetry [|13j as mq — > oo, and to quantify the importance of 
heavy-quark-symmetry-breaking pair creation effects residing in the 8E\ compared to 



their valence potential model counterparts [14 |. 

Finally, I note that while this paper in couched in the language of the nonrelativistic 
quark model, there is nothing in the proposed general framework that would prevent its 
being transferred to either a relativistic quark model or to field theory. 



13 



REFERENCES 



[1] P. Geiger and N. Isgur, Phys. Rev. D 41, 1595 (1990). 

[2] P. Geiger and N. Isgur, Phys. Rev. D 44, 799 (1991); Phys. Rev. Lett. 67, 1066 (1991); Phys. 
Rev. D 47, 5050 (1993); P. Geiger, ibid. 49, 6003 (1993); P. Geiger and N. Isgur, Phys. Rev. 
D55, 299 (1997). 

[3] Lowest order meson loop Feynman graphs can have two distinct topologies at the quark level: 
those with a closed qq loop and those where only the (relativistically) propagating valence lines 
appear. Both types of graphs, when analyzed in a hadronic basis, have mesons "in the air". In 
mesons, the latter class of graphs either correspond to rf meson loops (and so are dual to Z 
graphs) or to direct "s-channel" u> — 0-mixing-type OZI violations ||. In baryons, such graphs 
correspond to a variety of meson loops, but all are dual to Z-graphs. Given these dualities, one 
can imagine that non-closed-gg-loop diagrams are subsumed into the valence quark model. We 
focus here on the pure qq loop processes which are not otherwise taken into account in the 
quark model. 

[4] S. Okubo, Phys. Lett. 5, 165 (1963); Phys. Rev. D 16, 2336 (1977); G. Zweig, CERN Report 
No. 8419 TH 412, 1964 (unpublished); reprinted in Developments in the Quark Theory of 
Hadrons, edited by D. B. Lichtenberg and S. P. Rosen (Hadronic Press, Massachusetts, 1980); 
J. Iizuka, K. Okada, and O. Shito, Prog. Th. Phys. 35, 1061 (1966); J. Iizuka, Prog. Th. Phys. 
Suppl. 37-38, 21 (1966). 

[5] For a review of the adiabatic approximation in the context of the flux tube model (based on N. 
Isgur and J. Paton, Phys. Rev. D 31, 2910 (1985) and Refs. g]||) see N. Isgur in Proceedings 
of Few Body Problems in Physics, AIP Conf. Proc. 334, ed. F. Gross (A. I. P., New York, 1995), 
p. 3. For some early discussions of incorporating the nonperturbative effects of qq loops (which 
were all in the context of cc spectroscopy and decays), see J. Kogut and L. Susskind, Phys. 



14 



Rev. Lett. 34, 767 (1975); E. Eichten et al, Phys. Rev. Lett. 36, 500 (1976); and E.C. Poggio 
and H.J. Schnitzer, Phys. Rev. Lett. 41, 1344 (1978). These discussions are all different, and 
each introduces essential elements of the adiabatic approximation adopted here. 

[6] The considerations of this paper are trivially extended to mesons Q1Q2 with arbitrary masses 
mi and ni2- 

[7] R. Kokoski and N. Isgur, Phys. Rev. D 35, 907 (1987). 

[8] L. Micu, Nucl. Phys. BIO, 521 (1969); Jonathan L. Rosner, Phys. Rev. Lett. 13, 689 (1969); 
E.W. Colglazier and J.L. Rosner, Nucl. Phys. B27, 349 (1971); W.P Peterson and J.L. Rosner, 
Phys. Rev. D6, 820 (1972); Phys. Rev. D7, 747 (1973); A. Le Yaouanc, L. Oliver, O. Pene, 
and J.-C. Raynal, Phys. Rev. D 8, 2233 (1973); Phys. Lett. B71, 397 (1977); ibid. B72 , 57 
(1977); W. Roberts and B. Silvestre-Brac, Few Body Syst. 11, 171 (1992); P. Geiger and E.S. 
Swanson, Phys. Rev. D 50, 6855 (1994); Fl. Stancu and P. Stassart, Phys. Rev. D 38, 233 
(1988); 39, 343 (1989); 41, 916 (1990); 42, 1521 (1990); S. Capstick and W. Roberts, Phys. 
Rev. D 47, 1994, (1993); Phys. Rev. D 49, 4570, (1994). 

[9] One could instead expand Eq. ( |l~4| ) in €a Q +ej 3 mQ ' ) +g 2 /2/i Q/ g. This would lead to a Yukawa-like 
nonlocality in \p ' — p\ which, using a local approximation to this nonlocality and an analogue to 
Eq. (|3l|), would give an alternative definition of /±y^ proved (yy This version of /\v^ proved (r) 
would not contain the adiabatic limit as uiq — > 00, and so probably converges more slowly 
than Eq. (||). 

[10] If the pair creation occurs by creating a gap of classical length 2m^/6 in the flux tube, then 
V = m m ? m — h 7 — -7—^ — 77— . In other models the pair creation might be intrinsically nonlocal 

' mQ+rriq (mQ+m q )br r o j 

(for an example, see Ref. Q); in such cases the definition of /\y^P roved requires that one 
make a local approximation of the form of Eq. ([H|) to the model's actual matrix element 

( a (m Q ) p(m Q ) J R qq | QQ ( f ) ) _ 

[11] The explicit considerations given here for mesons are trivially extended to baryons and to 



15 



both meson and baryon hybrids. There is no reason to expect exact equality of the "improved" 
string tensions in these diverse systems. Indeed, S. Capstick and N. Isgur, Phys. Rev. D34, 
2809 (1986) note that the baryon string tension seems to differ by about 20% from the meson 
string tension. 

[12] See, for example, the unified description of mesons in S. Godfrey and N. Isgur, Phys. Rev. 
D32, 189 (1985). 

[13] N. Isgur and M.B. Wise, Phys. Rev. Lett. 66, 1130 (1991). 

[14] For an initial study of the relative importance of threshold and valence effects on spectroscopy, 
see N. Isgur, Phys. Rev. D 57, 4041 (1998) where the influence of low-lying S- wave thresholds 
on the sj e = ^ + and | + heavy quark mesons is considered. 



16 



